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Abstract: An analytic result for the 0(a s ) corrections to the triple differential 
(D \ B — > A u / i?j decay rate is presented, to leading order in the heavy-quark expansion. 

This is relevant for computing partially integrated decay distributions with arbitrary 
cuts on kinematic variables. Several double and single differential distributions are 
derived, most of which generalize known results. In particular, an analytic result 
for the 0(a s ) corrections to the hadronic invariant mass spectrum is presented. The 
effects of Fermi motion, which are important for the description of decay spectra close 
to infrared sensitive regions, are included. The behaviour of perturbation theory in 
the region of time-like momenta is also investigated. 
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1. Introduction 

Inclusive semileptonic decays of B mesons are of prime importance to determining 
the parameters \ V c b\ and \ V u b\ of the Cabibbo-Kobayashi-Maskawa matrix |T[. They 
also serve as probes for physics beyond the Standard Model, such as an extended 
Higgs sector [@] or right-handed weak couplings of the b quark [§, || . The total decay 
rates for these processes can be calculated in a systematic expansion in inverse powers 
of m h 0, |6], 0, ||, ||, 0, |ll[|. The same formalism can also be applied to calculate 
differential decay distributions, provided a sufficient sampling of hadronic final states 
is ensured by kinematics. Close to the boundary of phase space, the heavy-quark 
expansion must be generalized into a twist expansion to account for the effects of 
the "Fermi motion" of the b quark inside the B meson [|T^, [T3|] . 
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The leading contribution in the l/m& expansion is given by the free 6-quark de- 
cay into partons, calculated in perturbation theory as a power series in a s . Several 
authors have computed radiative corrections to various semileptonic decay rates and 
spectra. In particular, the QCD corrections to the total inclusive B — > X u I v\ decay 
rate have recently been computed to 0(a 2 ) [0, and the rate for B — > X c l&i is 
known to the same order from an extrapolation of exact results obtained at three 
different values of the invariant mass of the lepton-neutrino pair [115]. However, to 
our knowledge no results for the 0(a s ) corrections to the fully differential B — > Xlv\ 
decay distribution have been published so far. (An early investigation of these cor- 
rections was performed in ||16|| , where results are presented in complicated equations 
involving one-dimensional integrals.) Whereas this distribution by itself is not of 
direct phenomenological relevance (because it does not contain sufficient averaging 
over hadronic final states to be realistic), it is a necessary ingredient in the derivation 
of predictions for inclusive spectra with arbitrary cuts on kinematic variables. 

Here we present analytic results for the fully differential decay rate and several 
double and single differential distributions for B — > X u li>i decays. Throughout, we 
work to leading order in the heavy-quark expansion, omitting corrections of order 
(A/rrib) 2 , which have been discussed by previous authors (see, e.g., || 0, ||, |9|, |T0|, |TT]|). 



The QCD corrections are calculated including terms that do not contribute in the 
limit of vanishing lepton mass, so that our results allow treating the case of decays 
into t leptons. Semileptonic decays into final states containing a charm quark will 
be discussed elsewhere. A technical complication arises from the presence of soft and 
collinear singularities, which do not cancel in the fully differential decay distribution. 
These unphysical singularities appear because the calculation of inclusive rates is 
performed using external states containing free quarks and gluons. By virtue of global 
quark-hadron duality, the inclusive partonic rates are dual to the corresponding 
hadronic rates if a sufficient averaging over many final states is performed. We set 
m u = and regulate the infrared singularities by introducing a fictitious gluon mass 
A. The limit A — > can be taken at the end of the calculation and leads to singular 
distributions at p 2 = 0, where p is the total parton momentum. Inclusive spectra 
obtained by integration over a range in p 2 are infrared finite. 

2. Hadronic tensor at next-to-leading order in a s 

All strong-interaction dynamics relevant to inclusive semileptonic decays is encoded 
in the hadronic tensor 

W^(p,v) = --lmT^(p,v), (2.1) 

7T 

where 

TVfe „) = -i J <1> X e«™>- (*MIT{W,(0)}|f»M) (2.2) 



is the forward scattering amplitude given by the 5-meson matrix element of the time- 
ordered product of two weak currents, with J M = M7 M (1 — 75)6. For the calculation 
of QCD corrections it is convenient to choose the 6-quark velocity v (which can be 
taken to coincide with the velocity of the B meson) and the total parton momentum 
p as the two independent variables characterizing the hadronic tensor. The total 
momentum carried by the leptons is q = m^v — p. 

Because v 2 = 1, the two independent kinematic invariants are v ■ p and p 2 . The 
most general Lorentz-invariant decomposition of the hadronic tensor contains five 
invariant functions Wi(v ■ p,p 2 ), which we define as 

Wiivip, v) = Wi(v- p, p 2 ) [p^Vy + pvVp - g^v-p- it^p p a v J - 
- W 2 (v ■ p,p 2 ) + W 3 (v- p,p 2 ) V^Vy + 
+ Wi(v -p,p 2 ) (pftVy + p^^) + W 5 (v -p,p 2 )p^Pu- (2.3) 

At tree level W\ = 25{p 2 ) and W^i = 0. The five invariant functions W{ suffice 
to calculate arbitrary semileptonic decay distributions, including the case where the 
mass of the charged lepton is not neglected. In general, these distributions can be 
written in terms of three independent kinematic variables. One common choice of 
such variables is the charged- lepton energy Ei, the total lepton energy v ■ q (both 
defined in the 5-meson rest frame), and the invariant mass q 2 of the lepton pair. 
Here we choose a different set of variables, because the hadronic tensor is most 
conveniently calculated in terms of v -p and p 2 . Besides, experimentally the neutrino 
cannot be detected, whereas the total invariant mass and energy of the hadronic final 
state can be reconstructed directly. In terms of the parton variables, these quantities 
are given by 

s H = p 2 + 2A v ■ p + A 2 , E H = v ■ p + A , (2.4) 

where A = Mb—vh^. Experimental cuts on the region of low hadronic invariant mass 
or energy have been suggested as efficient ways to discriminate the small B — > X u I v\ 
signal against the much larger background from B — >• X c lvi decays |17], 18 , 0, [20], ^1 



Such a discrimination is important for a reliable determination of the Cabibbo- 
Kobayashi-Maskawa matrix element \V u b\. 

Although our results allow to treat the more general case, for the purpose of our 
phenomenological discussion we will set mi = 0. We introduce the scaling variables 

2£z -2 P 2 2v-p 

x = , p =— , z = , (2.5) 

nib m t m b 

in terms of which the triple differential decay rate is 

^^ = 12r {(l + x-.)(,-x-^)^^ 1 + (l-, + p 2 )^H/ 2 + 

+ [x(z -x)- p 2 } ^(W 3 + 2m b W 4 + m 2 b W^ , (2.6) 



where x — 1 — x, and 



r 



G 2 F \V ub \ 2 ml 
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The phase space for these variables is 

< x < 1 , x < z < 1 + x , max(0, z — 1) < p 2 < x(z — x) 



(2.7) 



:2.s) 



For fixed values of z andp 2 the lepton energy can vary in the range \{z— yj z 2 — Ap 2 ) < 
% < \{z + \/z' 2 — Ap 2 ), and since the hadronic tensor is independent of x it is possible 
to integrate over this variable to obtain the double differential spectrum 



1 dT 

r dzdp 2 



2 

2^2-4^1^(3 - 2z) - p 2 (A - 3z)) + 



+ ^l- z + p 2 ) r ^W 2 + 



+ (z 2 - Ap 2 [ 



(W 3 + 2m b W 4 + m 2 b W b 



where 



< z < 2 , max(0, z - 1) < p 2 < 



(2.9) 



(2.10) 





It is worth stressing at this 
point that any perturbative de- 
scription of inclusive decay rates 
must necessarily be formulated in 
terms of scheme-dependent para- 
meters such as the 6-quark mass 
m b and the related parameter A. 
These parameters are well-defi- 
ned to a given order in perturba- 
tion theory. In the context of our 
perturbative analysis, nib is to be 
identified with the one-loop pole 
mass of the b quark. The concept 
of a pole mass becomes ambigu- 
ous, however, if one goes beyond perturbation theory, and in that sense neither mj 
nor A are physical quantities. We will see in section ^ that, for normalized inclu- 
sive decay spectra, all reference to these parameters disappears once the leading 
nonperturbative corrections are included in the heavy-quark expansion. 

The 0(a s ) corrections to the hadronic tensor are obtained by evaluating the 
contributions of all physical cuts of the diagrams shown in figure [l], supplemented 
by wave-function renormalization graphs for the external b quarks. The sum of 
all contributions is gauge independent and free of ultraviolet divergences. There 
are, however, infrared divergences for p 2 — > from soft and collinear gluons, which 



Figure 1: One-loop diagrams contributing to the 
forward amplitude T fMU (p,v). 



we regularize with a fictitious gluon mass A = m^A. The limit A — > is taken 
whenever possible. We expand the invariant functions in a perturbation series as 
(with C F = 4/3) 



W i (z,i?)=W? > \z,?) + 



C F a s 

47T 



W?\z,f) + 0(ai), 



(2.11) 



where ^ = S(p 2 ), and all other functions vanish at leading order. Our results 
for the next-to-leading coefficients can be presented as follows: 

rn 
~2 



m " W [ l) = -5{p 2 ) 



41n 2 2 - 6 hi2 + + 4L 2 (1 - 2) + tt 2 + — 

1 — 2 2 



- /lR (A 2 ,2,p 2 ) + - 
8-z)(2 -2) 



i.ni±i + .nC 



+ 1 



2 t 2 



t 



+ 



l-t 



+ 



2 (B-z)(2-z) 



2z 



2 2 



4 3 



z 32 - 82 + 2 2 
~ + IzT 2 



2zH 2 
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+ ■ 



1, 1 + t 
- in , 

* l-t 



8zt 2 



1, 1 + t 
- In , 

t l-t 



8-32 32 + 22z-3z 2 3(12 - z) 
- — + — v — ; + 



Azt 2 



8t 4 



+ 



zt 2 5(4-2) 64 + 562- 72 2 3(12-2 



mi 
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16 16 

z In z 
+ 1 



+ 



Kf) - 1 - 



162t 2 16t 4 
32-52 3(12-2 



1, 1+t 

- In , 

t 1 -t ' 



2zt 2 



22t 4 



- 32 22 - 32 3(12 - 2) 
+ 



4 5 



42 22t 2 
1-22, 

In 2 - 



42t 4 



1, 1 + t 

- in , 

t 1 -t 



S(p 2 ) 



8 + 2 3(12-2) 



+ 



1 



42 22 2 t 2 



2 - 2 3(12 - 2) 



42 2 t 4 



22 2 t 2 

1, 1 + t 

- in 

t l-t 



22 2 t 4 



+ 



(2.12) 



Here t = — Ap 2 /z 2 , L 2 (x) is the dilogarithm, and we have omitted the step 
function 6{p 2 ) multiplying all regular terms. The expressions given above are regular 
for t — > 0, corresponding to the kinematic limit where p 2 — > 2 2 /4, i.e. |p | = 0. 

The function /ir(A 2 , z,p 2 ) entering the expression for w[ 1 ^ contains all reference 
to the infrared regulator. It is given by 

4(p 2 - A 2 ) 



f m (X 2 ,z,p 2 ) = 5(p 2 ) ln 2 A' + (5-41n2)lnA 2 + 9(p - A ) 



{p 2 - A 2 ) 2 + 2 2 A 2 



+ 



+ 



6{p 2 - A 2 

p2 



*!] [3 + ^] +41n[- + ^ 

P 2 / V P 2 J \P 2 Z 2 y 



(2.13) 
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In deriving this expression we have kept the regulator in all terms that diverge 
stronger than a logarithm in the limit where p 2 — > 0. The terms proportional to 
S(p 2 ) come from virtual corrections to the leading-order diagram, whereas the re- 
maining terms arise from cut diagrams with the emission of a real gluon. The limit 
A — > can be taken either if p 2 > by virtue of some experimental cut, or if the 
decay distribution is integrated over some range in p 2 . Studying the integrals of the 
function /ir with arbitrary regular test functions F(p 2 ), we find that in the sense of 
distributions one can replace /ir(A 2 , z,p 2 ) by 



f m (0,z,p 2 )=6(p 2 



f , 2 , 5\ A /hip 2 



(2.14) 



where the * distributions are defined as 

~6(p 2 - e) 









= lim 


\P 2 J 


* 


\np 2 ^\ 


= lim 


P 2 J 


* 



6(p 



+ 5{p 2 ) lne 



C ' ' s2 i 1 r/ -2\ i..'2 



hip 2 + -5(p 2 )hi 2 e 



(2.15) 



This definition is such that 



rrh 2 / 1 \ rrh 2 

^ dp 2 F{p 2 )\^j =F(0)lnm 2 + y o dp 2 



F(p 2 ) - F(0) 

p2 



(2.16) 



and similarly for the second distribution, so that the * can be omitted if the distri- 
butions are multiplied by a test function of 0{p 2 ). 

Using the result fl2.14| ), the 0(a s ) corrections to W\ can be rewritten as 



m t „,(i) 



Wl 1 



-KP 2 

-< 

+ i- 



8 ln 2 z - 10 In z + + 4L 2 (1 - z) + — 



1 



hip 2 



+ (81n*-7) - +- 



4 



1 - t 



+ 



-*)(2-*) 



+ 



2-z (8- z)(2 - z) 



1, l+t 
- in 

t 1 - t 



(2.17) 



2z ' 2zH 2 

The expressions (|2.12|) and (|2.17|) are the basis for all results presented in this paper. 



3. Double differential distributions 

The exact results for the invariant functions Wi given above allow the calculation 
of arbitrary B — ► X u lVi decay distributions to next-to-leading order in a s . In par- 
ticular, experimental cuts on the variables Ei, Eh and sjj can be implemented in a 
straightforward way if the distributions are obtained from a numerical integration 
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of (|2.6D . In this section, we derive analytic results for the double differential distri- 
butions in the variables [z,p 2 ) and (z,x) obtained after one integration of the fully 
differential decay distribution. These results are the basis for, e.g., the study of the 
hadronic energy spectrum with a cut on the charged- lepton energy (or vice versa), 
or the hadronic invariant mass distribution. The latter will be discussed in section [5| 

3.1 Distribution in the variables (z,p 2 ) 

Inserting the results for the invariant functions Wi into the general relation ( |2.9| ), we 
obtain a remarkably simple expression for the double differential decay rate. Defining 



i d 2 r 

r dzdp 2 



2z 2 (3-2z) 



6(p 2 ) + ^E 1 (z,f) 

All 



C F a s 
An 



E 2 (z,p 2 ), (3.1) 



where here and below we omit writing 0(a 2 ) for the neglected higher-order contri- 
butions, we find 



E 1 (z,f) = -S(p 2 ) 



8\n 2 z - 101nz + 



2\nz 



+ 4L 2 (l-z) + 5 



Air 2 



-4 



'lnp 2 \ 

E 2 (z,p 2 ) = ^P^5(p 2 



+ (81nz- 7) 



+ 



p A 



8 In— +7(1- t) 



2^(3 - Az) - 3(1 - 2z)p 2 - 2p 4 



, 1 + t 
In h 

1-t 



(3.2) 



1-z 

+ Azt (10 - 15z + 8p 2 ) . 
The kinematic range for the variables z and p 2 is given in (p . 1 0|) . 
3.2 Distribution in the variables (x, z) 

Another useful distribution is obtained by integrating the triple differential decay rate 
over the variable p 2 in the range specified in (|2.8|) . This leaves x and z as kinematic 
variables, which allows us to compute arbitrary distributions in the charged-lepton, 
neutrino or hadronic energies. The result for this distribution takes a different form 
for the two cases z < 1 and z > 1. For the first case, we find 



1 dT 



T dxdz 



12(2-x-z)(x + z- 1) 



Z<1 



1 + 



CpOL s E < (x,z 



An 2 — x — z 



(3.3) 



where 1 — z < x < 1, and 



F<(x,z) = -2(2- x - z) 



In' 



x + z — 1 



X 



+ 2L 2 (l-^) + 



2tt 2 



2(5 -2x- 2z) ln(l - x) - ^ In — 1 + 

15 1 x 



+ 



__ A 

15(x + z - 1) " 1 - x 30 



In 



(3.4) 



For the case z > 



1, we have instead 



1 d 2 r 

T dxdz 



Z>1 



C F a s 
Air 



F>{x,z, , 



(3.5) 



where < x < 2 — z and 



F>(x, z) = -24(2 - x - z)(x + z - 1) 



ln z 



x + z — 1 
1 — X 



- \n\z - 1) 



In fljy) and 

h 



p- (x + z — 1) ln(x + z — 1) — In ^ — - + 

5 5 1 — x 

+ ^(2-x-z) 3 \n(l-x) + ^ (2 - x - z) . 
, the coefficients fi are polynomials in x and z given by 



(3.6) 



71 - 44x + 26x 2 - 9x 3 + x A ) - (44 - 52x + 27a; 2 - 4a; 3 ) 2 + 
+ (26 - 27a; + 6a; 2 ) z 2 - (9 - 4x) z 3 + z 4 , 
/ 2 = QOx - 15 (2 - 3x + 5x 2 ) 2 + 5 (7 - 8x + 2a; 2 ) z 2 - 5 (2 - x) z 3 + z 4 , 
/ 3 = 15 (2 + 13a; - 7x 2 + 3a; 3 - a; 4 ) + (37 - 126a; - 29a; 2 - 32a; 3 ) z - 

- (1 - x)(:L7 - 7x)z 2 + 2(1 - x)z 3 , 

U = (71 - 115a; + 40a; 2 ) - 5 (23 - Ax) z + 5 (l4 - 9a; + 15a; 2 ) z 2 - 

- 5 (7 - 82 + 2x 2 ) z 3 + 5 (2 - x) z A - z 5 , 
f 5 = (l + Ax - x 2 ) + 2 (2 - x) z - z 2 , 

/ 6 = 15 (10 - 6a; + 2x 2 - x 3 ) - 2 (32 + 23a; + 16a; 2 ) z + 7 (2 - x) z 2 - 2z 3 . (3.7) 



In the limit x — > 1, we obtain the simple expression 



1 dT 



Tn dxdz 



12z(l - 2) 



x->l 



(3.i 



where 



Fi(x, z) = ln 2 (l - a;) - 21nzln(l - x) + - ln(l - x) + 



+ In z In z + 



In 2 
1-z 



+ 2L 2 (l-z) + 



2tt 2 



0(1 -x). (3.9) 



This result, which was previously derived by Akhoury and Rothstein ||22|| , is needed 
for the next-to-leading order resummation of Sudakov logarithms to all orders of 
perturbation theory pi . 



8 



4. Single differential spectra 



While the results presented in the previous section are new, there already exist several 
calculations of the 0(a s ) corrections to single differential spectra in B — > X u lv\ 
decays. Here we derive the distributions in the variables x, z and p 2 . This allows 
comparison of our results with existing calculations. 

4.1 Charged-lepton energy spectrum 



Integrating the double differential decay rate derived in section |3.2j over z yields the 
spectrum in the variable x, which measures the energy of the charged lepton in the 
I?-meson rest frame. We obtain 

C F a s 



r ax 



2x) 



1 - 



2tt 



G{x) 



(4.1) 



where < x < 1, and 
G{x) 



in 2 n 



x 



j- z \ 2tt 2 82 - 153s + 86x 2 

+ 2L 2 (x) + — + — + 



+ 



41 - 36x + 42x 2 



3 

16a; 3 



12x(3 - 2x) 



ln(l 



x 



6x 2 (3-2x) 

This agrees with the well-known result obtained first by Jezabek and Kiihn 

The function G(x) is regular for x— > 
0. On the other hand for x— >1, i.e. close 
to the boundary of phase space, there are 
Sudakov logarithms reflecting the incom- 
plete cancellation of infrared divergences 
due to the limited phase space available 
for real gluon emission: 



(4.2) 



51. 



31 

G(x) = ln 2 (l -x) + — ln(l - x) + 



7T 



6 

; + 0(l 




X) 



(4.3) 



These endpoint singularities are integra- 
ble, and the total decay rate is given by 

Cpot 



Figure 2: Lepton energy spectrum dr/dx 
in units of T, at tree level (dashed) and in- 
cluding 0(a s ) corrections (solid). 



2n 



7l~ 



25\ 

t) 



(4.4) 



In figure |2|, we show the result for the charged-lepton energy spectrum obtained at 
leading and next-to-leading order in perturbation theory, using a s = 0.22. Here and 
below we normalize the distributions to the total decay rate T, so that the spectra 
shown have unit area. At tree level we use T = r , whereas at next-to-leading order 
we take the result given in ( f4.4| ). It is evident from the figure that the perturbative 
corrections affect the spectral shape close to the endpoint only. 
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4.2 Hadronic energy spectrum 

Because of the relation 
Eh — A 



v ■ p 



Mb-W+E^-A, (4.5) 



both the hadronic energy spectrum and 
the distribution of the total energy of the 
lepton pair can be derived from the dis- 
tribution for the scaling variable z, which 
is obtained by integrating the double dif- 
ferential distribution in (|3.1|) over p 2 . Be- 
cause of the phase-space constraint 
shown in ( [2.10| ) , the resulting expressions 
are different for the two cases z < 1 and 
z > 1. In the second case p 2 = is not 
allowed by kinematics, and thus only di- 
agrams with real gluon emission contribute. We find 




Figure 3: Hadronic energy spectrum dr/dz 
in units of V, at tree level (dashed) and in- 
cluding 0(a s ) corrections (solid). 



i dr 

r dz 
i_ar 

r dz 



2z 2 {3 

CfCVs 

2tt 



- 2z) 
H>(z), 



1 



CfCXs 

2tt 



< z < 1 



1 < z < 2 



(4.6) 



with 



2L 2 (1 - z) + tt 2 + 



2z 2 (3-2z) 



ln 2 (z 



9-Az 
3-2z 

- 1) - 



In z 



4860 - 3720z + 585z 2 - 42z 3 + Az 4 



2\n 2 z 





" 2 1 








+ 



(2 - z) (1248 + 3798,2 - 2946z 2 + 5172 s 



360(3 - 2z) 
+ 

34z 4 + 4z 5 ) 



+ 



12z + 12z 2 



180 
ln(z- 1). 



+ 



(4.7) 



Our results for the functions H K (z) and H > (z) agree with the findings of Czarnecki, 
Jezabek and Kiihn [p6 |. Whereas the function H < (z) is regular for z — > 1, H > (z) 
exhibits a logarithmic divergence at this point, because the singularities from soft 
gluon emission are not compensated by virtual gluon corrections. We find 



H>(z) =21n 2 (2-l) + 71n(z-l) 



2tt 2 2587 ^. 

+ ^rr + 0(z - 1) . (4. 



3 180 

The result for the hadronic energy spectrum at leading and next-to-leading order in 
perturbation theory is shown in figure ^. The double-logarithmic singularity at the 
point z = 1 located inside the allowed kinematic region provides an example of a 
"Sudakov shoulder" J27 |. 
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4.3 Parton mass spectrum 



Because the parton invariant mass y/p 1 is not an observable quantity, the spectrum 
in this variable is not of direct phenomenological relevance. However, in the region 
p 2 = 0(m 2 ) it follows from ( |2.4| ) that p 2 ~ sh up to corrections of order A/Mb, and 
thus the spectrum dT/dp 2 is a reasonable approximation to the hadronic invariant 
mass distribution. Also, this spectrum can be used to compute moments of the 
hadronic mass distribution, since for instance (s#) = m 2 , (p 2 ) + Amf 1 (z) + A 2 . Finally, 
as we will discuss in section |7|, the parton mass spectrum allows us to perform a 
quantitative study of the behaviour of perturbative QCD in the region of time-like 
momenta. 

Integrating the double differential decay rate in ( |2.9| ) over z in the range 2^/p 1 < 
z < 1 + p 2 , we obtain 



i dr 2 

7T = $(P ) 

r dp 2 



i 



C F a s 
2tt 



7T~ 



187 
~T2 



2vr 



I{p 2 



(4.9) 



where < p 2 < 1, and 




6p 2 + ^p 6 ) \nf 



(4.10) 



The regular terms in this result, i.e. the terms that remain for p 2 > 0, have been 



computed previously by Falk, Luke and Savage [28|, and we agree with their result. 



5. Hadronic invariant mass distribution 

Imposing a kinematic cut sh < Mp on the inclusive semileptonic decay rate of B 
mesons is an efficient way to separate the Cabibbo-suppressed signal from b — » u 
transitions from the background of b — > c decays [|T7]]. Using the relation between 
the parton variables p 2 and v ■ p and the hadronic invariant mass displayed in ( |2.4j ) , 
and denoting §h = sh/tu 2 and e = A/m&, we find that 

dr d 2 r 

dsn Jz Q dz dp 2 

where 

z = 2(\fhi - e) , zi = min^ Sg g £ , l-e + ' ( 5 ' 2 ) 

Because of the form of z\, one must distinguish the cases where sh is smaller or larger 
than e(l + e). In the second case, only diagrams with real gluon emission contribute 



e 2 < s H < (l+e)' 



(5.1) 



n 



Sh 



EZ 
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to the spectrum. We define 



' <lJ T(s H ,e)- < =^J < (s H ,e), e 2 < s H < e(l + e) , 



T ds H 
1 dr C F a 



2tt 

J>(s H ,e) 



e(l + e) < s H < (1 + e) 



T dSff 27T 

The exact results for the functions entering this expression read 
eT{s H ,e) = 2x 2 {3 - 2x) , 



(5.3) 



£ J<(sr, s) = 2x 2 (3 — 2x) 



vr 2 + 2L 2 (l-x) -2L 2 (-^) 



1 f x 

+ q — ^2 — 2x 2 d 3 )(x + e) In 1 1 H — 



+ 
+ 



+ 2a; 2 (9 - 4e) In a; - - edi + - rf 4 - — d, 



x 



X 



£ J>(s H ,e) = 2x 2 (3 - 2x) x 

J2 



X 



7T 



— + ln z y - 2 hTz + 4L 2 — 



4Lo - -2L 



+ - {sd\ — x<i 2 — 2x 2 d 3 )(x + e) In 
+ (7 + Ay - V — d e + ^ d 7 )\ny + 



y_ + e_ 

l + e 



+ 



3e 3 

I — y 

+ TTTT. — r(^8 + %dg — 2x 2 d\o) 



18(1 + e) 



where 



x 



s H - e< 



y 



s H -e(l + e) 



(5.4) 



(5.5) 



: 1 + £ 

and di are polynomials in e given in equation ( |A.l ) of the appendix. 

Of key importance to the determination of \V u b\ is the question which fraction 
of all B — > X u I V\ events have hadronic invariant mass below the charm threshold. 
To address this issue, we compute the integral of the spectrum up to a cutoff m 2 = 

e 2 < m 2 < (1 + e) 2 , (5.6) 



T(m 2 , e) 
and write the result in the form 
T(m 2 ,e) = T 



dr 

ds H 



ds H 



./ a 2 \ CpOl s ■ / - 2 \ 

H m , £ ) o — -?<( m ' £ ) 

Z7T 



r(m 2 ,e) = r c 



C F a s . . „ 2 
1 = — J>{m ,ej 



Introducing the variables 



2tt 



m 2 — e 2 



e 2 < m 2 < e{l + e) , 
e{l + e) <m 2 < (l+e) 2 . (5.7) 



rh 2 — e(l + e) 
l + e 



(5.f 
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0.25 




0.25 



m 



Figure 4: Left: Invariant hadronic mass spectrum dT/dsn m units of T, at tree level 
(dashed) and including 0(a s ) corrections (solid). Right: Partially integrated spectrum 
T(rh 2 ,e) in units of T. 



we obtain 

t{m 2 ,e) ■ 
j<(m 2 ,e) : 



/i 3 (2-/i), 
/i 3 (2-/i) 

■2 



7T~ 



2L 2 (l-/i)-2L 2 - 



+ 



7r" / 20 5 \ 

+ — - 2L 2 (1 - fi) + ^2// + /i 2 + — /i 3 --// 4 J ln/i + 

+ ^ (^ 2 ei + /iee 2 - /i 2 e 3 + /i 3 e 4 ) (// + e) ln^l + ^ 

/i ji 2 /i 3 /i 4 
-6 e5 -12 e6 -18 e7+ 72 e8 ' 



-/i 3 (2-/i) 
f 7r 2 + ln 2 g 



y + ln 2 ^ - 2 ln 2 /i + 4L 2 (-i) - 4L 2 (- 



+ 



1 



6e(l + e) 



6 



e 9 + /ie 10 + /i 2 eeii - /i 3 e i2 ) gin g 
g + £ 



— — [e ei + /i ee 2 — n e 3 + [i e 4 ) (/i + e) In 



1 + e 



+ 



1 + e 



72 



ei3 



e 14 



/i 2 /i 3 /i 4 e N 
12 615 + 18 616 "72" 617 , 



(5.9) 



with coefficients given in equation ( |A.2j ) of the appendix. 

In figure we show the results for the hadronic invariant mass spectrum (left- 
hand plot) and for the integral of this spectrum up to a cutoff m 2 (right-hand plot) 
for e = 0.1, corresponding to nib ~ 4.8 GeV. Only the lower portion of the kine- 
matic range for the variables sh and rh 2 is displayed. In contrast with the energy 
spectra considered in the previous section, radiative corrections have an important 
impact on the shape of the hadronic invariant mass spectrum and lead to a signif- 
icant redistribution from lower to higher masses. Nevertheless, it is apparent from 



o 



O 
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the right-hand plot that imposing a cut sh < Mp, corresponding to m 2 « 0.15, 
would leave most of the B — > X u Z events unaffected but at the same time remove 
all B 



X c I v\ events [1 7, 20 



] . Because such a cut falls close to the sharp edge of 
the perturbative hadronic mass spectrum, however, a careful treatment of nonper- 
turbative corrections is necessary before any conclusions can be drawn. This will be 
discussed in the following section. 



6. Implementation of Fermi motion 

The perturbative results presented above are only one ingredient to a consistent the- 
oretical description of inclusive decay spectra. In kinematic regions close to phase- 
space boundaries these spectra are infrared sensitive and receive large nonpertur- 
bative corrections. Because the corresponding effects can be associated with the 
motion of the b quark inside the B meson, they are commonly referred to as "Fermi 
motion". These effects are always important when the perturbative prediction for 
an inclusive decay distribution exhibits a rapid variation on a scale that is paramet- 
rically smaller than nib- Usually, in such a case one encounters large perturbative 
corrections from Sudakov logarithms. For the single differential spectra considered 
in this paper, this happens in the endpoint region 1 — x — 0(A/ra&) of the charged- 
lepton energy spectrum shown in figure [|, the central region 1 — z — 0(A/mb) of the 
hadronic energy spectrum shown in figure ^|, and the low- mass region §h = 0(A/m&) 
of the hadronic invariant mass spectrum shown in figure [|. Note that in the first 
two cases these are small fractions of the kinematic regions; however, in the case of 
the hadronic mass distribution essentially all of the spectrum is concentrated in the 
region §h = 0(A/m&), where nonperturbative effects are important. 

Fermi motion effects are included in the heavy-quark expansion by resumming 
an infinite set of leading- twist corrections into a shape function F(k + ), which governs 



the light-cone momentum distribution of the heavy quark inside the B meson [12 



T3j. The physical decay distributions are obtained from a convolution of parton 
model spectra with this function. In the process, phase-space boundaries defined by 
parton kinematics are transformed into the proper physical boundaries determined 
by hadron kinematics. The shape function is a universal characteristic of the B 
meson governing inclusive decay spectra in processes with massless partons in the 
final state, such as B — > X u I V[ and B — > X^. The convolution of parton spectra 
with this function is such that in the perturbative formulae for the decay distributions 
the 6-quark mass mb is replaced by the momentum dependent mass nib + k+, and 
similarly the parameter A = Mb — rrib is replaced by A — k + |fl2|] . Here k + can take 
values between — and A, with a distribution centered around k + = and with a 
characteristic width of 0(A). Introducing the new variable q + = A — k + , it follows 
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that, e.g., the scaling variables x and z are replaced by the new variables 

x = 2El z = 2 ^ ~ g±j (61) 

q Mb — q+ Zq M B -q + ' ^ 

and the physical spectra for the charged-lepton energy and for the total hadronic 
energy are, respectively, given by 

dr f M B -2 El F(A~q + ) dr M B 
di?; Jo Mb — g + dx 

and 



dr F(A - o+) dr 

- 1 ' dg + 4 ^-(z q ) , < E H < M B . (6.3) 



dE H Jo M B — q + dz 

The perturbative spectra dT/dx and dT/dz have been given in ( |4.1| ) and ( |4.6| ). The 
upper limits of the q + integration follow from the allowed kinematic ranges for the 
variables x q and z q . Similarly, for the hadronic invariant mass spectrum we define 



Sh « - {Mb ~ <1 + Y ' _ M B — q+ (6 ' 4) 



and obtain 



dr rv^F F(A-q+) dr /a2 n ,, 2 , . 

j = / dg + -A ^±1 s 2 £g , < s H < M 2 B , 6.5 

ds H Jo (M B - q + Y ds H 



with dr /dsu as given in ( |5.3|) . Here the upper limit of the q + integration is enforced 
by the requirement that s Hq > e 2 q . rrom ( |6.5| ), it follows that the integral over the 
hadronic invariant mass spectrum from to some cutoff s^ ax takes the form 



T(^ ax ) = ds H -5L = / V dq + F(A - g + ) r(m , e„) , (6.6) 
Jo ds H Jo H 

where m 2 q = s H ax /(M B — q+) 2 , and the quantity T(m 2 ,e) has been defined in ( |5.6|) . 

After the implementation of rermi motion the kinematic variables take values 
in the entire phase space determined by hadron kinematics, rather than the parton 
phase space appropriate for perturbative calculations. R)r instance, the maximum 
lepton energy attainable is Ef 13 ^ = Mb/2 rather than m&/2. In the above formulae 
the only reference to the nonperturbative parameter A resides in the shape function. 
All other mass parameters refer to physical hadron masses or energies. We stress 
that for the numerical evaluation of the convolution integrals it is necessary to have 
explicit analytic results for the perturbative spectra. R)r the important case of the 
hadronic invariant mass distribution, such results have not been presented so far in 
the literature. 1 



1 In |2l| one can find a figure showing numerical results for the function r(m 2 ,e) for some 
particular choices of parameters. These results do not permit the computation of the physical 
quantity T(s™ ax ). 
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Figure 5: Left: Charged-lepton energy spectrum dT/dEi (in units of T x GeV - ) versus 
Ei for nit, = 4.8 GeV (solid), 4.65 GeV (red dashed), and 4.95 GeV (green dashed). Right: 
Same for the hadronic energy spectrum dT/dEn- 



Several functional forms for the shape function have been suggested in the liter- 
ature. They are subject to constraints on the moments of this function, A n = 
which are related to the forward matrix elements of local operators on the light 
12fl . The first three moments satisfy Aq — 1, A± — and A 2 = where f/£ 



cone 



is the average momentum squared of the b quark inside the B meson [29|. For our 
purposes, it is sufficient to adopt the simple form pO 



F(k+) = N (1 - x) a e {l+a)x ■ x=i±<l, (6.7) 

which is such that A\ = by construction (neglecting exponentially small terms in 
mj/A), whereas the condition A = 1 fixes the normalization N. The parameter a 
can be related to the second moment, yielding A 2 = |/i 2 = A 2 /(l + a). Thus, the b- 
quark mass (or A) and the quantity a (or ^ ) are the two parameters of the function. 
A typical choice of values is m;, = 4.8 GeV and a = 1.29, corresponding to A w 
0.48 GeV and fi^ ~ 0.3 GeV 2 . Below, we keep a fixed and consider the three choices 
mb = 4.65, 4.8 and 4.95 GeV. The spread of the results provides a realistic estimate 
of the theoretical uncertainty associated with the treatment of Fermi motion. This 
uncertainty could be removed if the shape function were extracted, e.g., from a precise 
measurement of the photon energy spectrum in B — > X s ^y decays |30| . 

Figure [5] shows the results for the charged-lepton energy spectrum and for the 
hadronic energy spectrum, including Fermi motion effects. The three curves in each 
plot correspond to different values of the 6-quark mass. Comparing the shape of 
the spectra with the perturbative results shown in figures ^| and |3| indicates that 
for the charged-lepton energy spectrum nonperturbative effects are important in the 
region Ei > 1.8 GeV, whereas they affect the hadronic energy spectrum in the range 
2.6 GeV < Eh < 3.6 GeV. Generally, Fermi motion effects smooth out any sharp 
structures in the perturbative spectra. From the results for the hadronic energy 
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Figure 6: Left: Hadronic invariant mass spectrum dT/dsH (in units of T x GeV -2 ) versus 
sh- The meaning of the curves is the same as in figure [|. In addition, the short-dashed 
line shows the spectrum for mj, = 4.8 GeV without Fermi motion effects included. Right: 
Fraction of B — > X u I i>i events with hadronic invariant mass below a cutoff s^ ax . 



spectrum it follows that the fraction of events with E H < M D is about (30 ± oj7o, 
with a moderate uncertainty from the dependence on m?,. Thus, imposing a cut on the 
energy of the hadronic final state in B — > X I v\ decays provides a reasonably efficient 
way of separating the rare b — > u transitions from the much larger background of B 
decays into charmed particles |TJ, [T9j . 

In Figure || we show the results for the hadronic invariant mass spectrum and 
for the fraction of B — > X u 1 1>\ events with hadronic mass below a cutoff s^ ax . In each 
plot, the short-dashed line shows for comparison the perturbative spectrum obtained 
with rrib = 4.8 GeV and ignoring the effects of Fermi motion. The difference between 
this curve and the solid one is only due to nonperturbative effects. Clearly, these 
effects have a very important impact on the shape of the spectrum in the entire 
low-mass region relevant to experiment, and even in the region close to or above the 
charm threshold. Another important observation is that Fermi motion effects remove 
completely the sharp structure at sh ~ A Mb (corresponding to §h ~ £(1 + e)) 
in the perturbative hadronic mass spectrum and replace it by a broad bump at a 
significantly lower value of s#. In other words, the value of the unphysical quantity 
A is of no direct relevance to the hadronic invariant mass spectrum. 

From the right-hand plot in figure || we deduce that the fraction of events with 
hadronic invariant mass below the charm threshold (s# < 3.49 GeV 2 ) is about (80 ± 
10)%. Therefore, applying a cut s# < M|, would be a most efficient discriminator 
between semileptonc b — ► u and b — ► c transitions, which would allow for a largely 
model- independent determination of \V u b\. If for experimental reasons the cutoff on 
the hadronic mass is lowered to (1.5 GeV) 2 , the fraction of contained events drops to 
about (60 ± 15)%, which is still significant. Even in this more pessimistic scenario, 
\Vub\ could be extracted with a theoretical uncertainty of about 15%. 
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7. Perturbative QCD in the time-like region 



As a last application, we use our results for the perturbative corrections to the in- 
clusive B — > X u I v\ decay distributions to investigate the behaviour of perturbative 
QCD in the region of time-like momenta. The heavy-quark expansion used in the 
calculation of inclusive decay rates relies on an application of the operator prod- 
uct expansion (OPE) in the Minkowskian region. In this region there are physical 
singularities from on-shell intermediate hadron states, which are not reproduced by 
perturbation theory. The hypothesis of quark-hadron duality is the assumption that 
the OPE can still be applied provided there is sufficient averaging over hadronic final 
states, so that the properties of individual hadron resonances become unimportant. 

There is at present no known way of quantifying from first principles how well 
this assumption holds in the case of inclusive heavy-quark decays. In this section, 
we investigate the question whether perturbation theory itself signals the problem 
by exhibiting singularities when the resonance region is approached. To this end, 
we express the differential decay rate in (z,p 2 ) as the imaginary part of a correlator 
T(z,p 2 ), 

^-^-~ 2 = -ImT(z,p 2 + ze), (7.1) 

1 o CLZ djr 7T 

where z is taken to be real and inside the interval < z < 2. Up to irrelevant 
numerical factors, this correlator is the contraction of the tensor T^iv^p) defined 
in ( |2.2j ) with the lepton tensor, integrated over x. The imaginary part is nonzero 
if p 2 is in the interval given in ( 2.10Q . It follows that the correlator satisfies the 
dispersion relation 

T{z,p) = - ds , (7.2) 

7T Jzq S — p Z 

where zq = max(0, z — 1). In a similar way, we can express the differential spectrum 
in p 2 as the imaginary part of a correlator T(p 2 ), 

-r^ = -ImT(p +ie) , T(p ) = - ds , 7.3 

1 o dp Z 71 7T JO S — p 

where the imaginary part is nonzero if < p 2 < 1. With this definition, it follows 
that 

T(p 2 )= f 2 dzT(z,p 2 ). (7.4) 

The correlators T(z,p 2 ) and T(p 2 ) as functions of complex p 2 contain information 
about the behaviour of perturbation theory close to the region of physical singu- 
larities. For the purpose of illustration, we will discuss the structure of T(p 2 ) in 
detail. An analoguous discussion could be made for the correlator T(z,p 2 ) at any 
fixed value of z. 
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Performing the dispersion integral in (7.3) using the results for the 0(a s ) cor- 
rections given in ( f4.9|) and (|4.10| ), we obtain 



where 



K{p 2 



2 + T p 



C F a s 
2tt 



K{p 2 



(7.5) 



6 18 ^ 18^ 18 y 1 1 



31 41 
" 18 
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-p 
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+ 



233 



79 



7T 



72 



p H p" 

36 F 18 F 



(7.6) 



The function K(p 2 ) is analytic in the cut p 2 plane with a discontinuity along the 
interval < p 2 < 1, as required by the analyticity properties of the correlator 
following from (\l.3\). Note that in perturbation theory the only singular point is 



p 2 = 0, where the function has a logarithmic singularity. To get a reliable answer for 
a physical quantity, the region close to this singularity must be avoided. 

Consider, as an example, the inclusive B — > X u I V[ decay rate integrated over an 
interval in p 2 including the origin: 



Tim 1 



dr 



2 — - - u 

o dp z 2m Jo 



dp 2 T(p 2 + ie) — T{p 2 — ie) 



(7.7) 



Using Cauchy's relation, the integration contour can be deformed into a circle in the 
complex momentum plane touching the real axis at the point p 2 = m 2 . We obtain 



Ylrhr 



r 



i - 



C F a s 



2tt Jo 2tt 



27r d( P 2 

— K[m e^) 



(7- 



The situation is illustrated in figure 0. Note that for m 2 > 1 the result of the 
integration is independent of rh 2 , and it is thus possible to take the radius of the circle 
arbitrarily large, so that the contour is far away from the singularities. Therefore, is it 
generally accepted that the total inclusive semileptonic decay rate can be calculated 
reliably in perturbation theory. 

If m 2 < 1, on the other hand, the contour probes the region of physical resonances 
at the point where it touches the cut. It is natural to ask whether perturbation theory 
exhibits singularities in the vicinity of the cut, which could signal a breakdown of 
quark-hadron duality in the calculation of partially integrated decay rates such as 
r(m 2 ). Remarkably, we find that this is not the case. As long as rh 2 is not too close 
to the origin, the perturbative corrections to the correlator T(p 2 ) are well behaved 
everywhere in the complex momentum plane, even close to the cut. This is evident 
from figure which shows contours of the real part of K(p 2 ) in the complex p 2 
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Figure 7: Deformation of the integra- 
tion contour in the evaluation of the 
quantity T(m 2 ). 



Rep' 



Figure 8: Contours of KeK(p 2 ) in the 
complex p 2 plane. The contour lines re- 
fer to values between 2.5 (innermost right 
contour) and 4.5 (innermost left contour) 
in units of 0.25. 



plane. 2 We take this as an indication that partially integrated decay rates can be 
calculated using the heavy-quark expansion as long as they are not restricted to a 
range in p 2 too close to the origin. 

8. Conclusions 

We have presented analytic results for the next-to-leading order perturbative cor- 
rections to the triple differential B —>■ X u I v\ decay rate. They provide the basis 
for the computation of arbitrary inclusive decay distributions to 0(a s ), including 
experimental cuts on various kinematic variables. Our results are sufficiently general 
to allow treating the case where the mass of the charged lepton cannot be neglected. 
As an application, we have presented explicit results for several double and single dif- 
ferential distributions, most of which had not been derived previously. In particular, 
we have discussed in detail the 0(a s ) corrections to the hadronic invariant mass dis- 
tribution, which are an important ingredient in a theoretically clean determination 
of the element \V u b\ of the quark mixing matrix. 

We have shown how the leading nonperturbative corrections affecting inclusive 
decay spectra can be incorporated in a QCD-based framework by convoluting the 
perturbative spectra with a 6-quark momentum distribution function. This is im- 
portant for addressing the question of how experimentally one may separate the 
B — > X u I i>i signal from the large background of semileptonic decays into charmed 



o 



O 



2 The imaginary part does not contribute to the integral along the circle in 
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particles. We find that (30 ±5)% of all B — > X u lv\ events have hadronic energy 
below the charm threshold, while (80 ±10)% have hadronic invariant mass below 
M|,. If the cutoff on the hadronic mass is lowered to (1.5 GeV) 2 , this fraction drops 
to (60 ±15)%, which would still allow for a largely model-independent determination 
of \V uh \. 

Finally, we have studied the behaviour of perturbative QCD in the complex 
momentum plane, finding that there is no evidence for large corrections except for 
the region close to p 2 = 0. This observation can be taken as circumstantial evidence in 
support of global quark-hadron duality, which underlies the heavy-quark expansion 
for inclusive decay rates. 
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A. Coefficients entering the hadronic mass spectrum 

Here we list the coefficients di and entering the exact theoretical expressions for the 
hadronic invariant mass spectrum in (|5.4|) , and for the integrated spectrum in ( |5.9| ). 
These coefficients are polynomials in e = A/rrib given by 





= 9 ± 10e - 7e 2 - 2Ae 3 - 40e 4 , 




d-2 


= 9 - 20e - 7e 2 - 6e 3 + 32e 4 , 




d 3 


= 1 + Se - 3e 2 + 2e 3 , 




di 


= 153 + 50e + 7e 2 - I2e 3 - 104e 4 , 




d 5 


= 167 + 34e - 39e 2 + 40e 3 , 




d 6 


= 3 + Se + Ae 2 , 




d 7 


= 1 - e - 6e 2 - 4e 3 , 




d s 


= 143 ± 172e + 113e 2 + 42e 3 - U6e 4 - 


- 384e 5 - 240e 6 


d 9 


= 348 + 497^ + l7Ae 2 - 7e 3 - 388e A - 


312e 5 , 


dio 


= 138 + 159e - As 2 + I7e 3 + 40e 4 , 





and 

d = 4 - I2e 2 - 25e 3 - 30e 4 , 

e 2 = 14 + 20e - 2e 2 - 23e 3 - 50e 4 , 

e 3 = 14 - lOe - 2e 2 - he 3 + 22e 4 , 
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